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o , 

^ By introducing the intrinsic pair operators which commute with number fluctuation 

^SJ , operator, a new formalism is given for the number-conserving description of the pairing cor- 

■ relations. The difficulty in the conventional RPA treatment for pairing rotation is removed. 

o 

■ It is well known that the pairing correlations play an important role in various 
I phenomena in nuclei The correlations are usually described by the BCS theory 

or Hartree-Fock-Bogoliubov (HFB) theory. In a super phase the particle number is 
conserved only in average and its fluctuation is important in a finite system such as 
nuclei . 

The projection method is the rigorous one to restore the particle number con- 
servation However it loses the simple picture of excitation modes. On the other 
hand the RPA method treats the number fluctuation mode as a zero-energy eigen 
mode and recovers the broken symmetry within its approximation . However it 
, has some difficulties related with the small amplitude approximation. For example, 

' the norm of the RPA ground state becomes divergent. Extension of RPA has been 

made in various ways although such an approach becomes rather complicated . 
O , In this letter we propose a new method to treat the number conservation. We 

^ I first define the intrinsic pair modes which commute with the particle number operator 

and extract the intrinsic Hamiltonian made up of only these modes. The wave 
function is formally given by the number projection on the intrinsic sates. However 
^ ' the integral calculation for the projection is not necessary, because the eigen value 

■ problem of the original Hamiltonian reduce to that of the intrinsic Hamiltonian. 
Thus we take account of both merits of RPA and the projection method. 

We start with a general effective Hamiltonian, 

H = ^ eacl^Ca + ^ ValSjSclc^/sCSCy, (1) 

where cj^ and Cq, are the creation and annihilation operator for the single particle 
state with the quantum number a and ea denotes the single particle energy. For 
simplicity we discuss the system with one kind of nucleon (proton or neutron) and 
the BCS case in the following. 

After the Bogoliubov-Valatin transformation, 

al = Uacl - VaCa, (2) 
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the particle number operator N = J2a (it^a is written by the quasi-particle operators 
as, 

N={N) + Nc + N„ (3) 
^c = Yl '^aVaia^al, + aada), = X](^a " vl)al^aa, (4) 

a a 

where a denotes the time reversal of a. {N) represents the expectation value of N 
for the BCS vacuum |0). The number fluctuation operator is defined as 

AN = N-{N)=Nc + Ns. (5) 

Nc is rewritten as 

N, = LiXl+X,), L' = {AN^) = 2Y,ulvl (6) 

a 

by introducing the Tamm-Dankoff (TD) type collective operator Xl, 

+ X - «««« 

Xl = Y,Mc^(i)^, (7) 



V;,(a/3) = EV'cM) = l- (8) 

^ a/3 

In super phase L is assumed to be large so that we regard £ = 1/L as a small 
parameter. Noticing that (0| XcXl |0) = 1, Xc and are regarded as a quantity 
of order and therefore Nc is 0(e~^). On the other hand 

(0| X.NsXl |0) = Y.{ul - vl)i^l{aa) = 1 -2Y,vl^l{aa) < 1 (9) 

a a 

SO that Ng is at most 0{e^). 

Let us define the TD type non-collective mode as 

44 



Xt=^^4a/3)^, (10) 



which satisfies 

(0|X,Xt|0) = 0, (0|X„X^,|0) = 5nn'. (11) 
This leads to the following equation for tpn{ctP) 

^ Vc(a/?)V'n(a/?) = 0, ^Mo^P)4^n'{aP) = Sun'. (12) 
a/3 aP 

Such a jpniciP) can be obtained, for example, from the eigen equation for stability 
matrix, 

^^(•^"PitS + ^af3;'yS)^n{lS) = Un^niaP) (^n > 0), (13) 
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where 

Aa0.,^5 = {O\[apaa,[H,a\al]]\O), Ba^.^s = {0\[[H,al,al],ala\]\0) (14) 

It is well known the eigenvector with zero eigenvalue corresponds to ipcioiP). It should 
be mentioned that the wave functions V-„(«/^) thus determined are not unique and 
there remains the freedom of unitary transformation between them. However they 
span the complete set together with X}. in the sense that 



"^ = EV'mM)4' (15) 



where we use the notation /j, = c oi n. 

Now we define the intrinsic pair operators Xj^ and Bq corresponding to and 
Bq = al^ap by the following equations; 

[AN,Xl]=0 Xn\0) = 0, (16) 

[AN,Bg] = 0, B,\0) = 0. (17) 
Here we assume that they are obtained by an expansion form, 

Xl = xl + exn'^ + e'xn'^ + ..., (18) 

B, = Bf +eSW+£25f + .... (19) 

The leading term Bg^^ is given by 

= Bq- r^{q)XlX^ - [rcn{q)XiX^ + r„e(g)Xtx„] , (20) 

n 

where r^^^q) is defined as 

r^,{q) = 2Y,^^{a-i)Mh)- (21) 

7 

Note that [Xc, Bq] = and [Xl,Bq] = in lowest order. Prom Eqs.(18) and (19) the 
pair operators X^ and Bq are conversely expanded as 

Xl = Xi + eXt + e^xt (2) + . . . , (22) 
S, = Sf + + ..., (23) 

where X\^'^^ {i > 1), Bq^ {i > 0) are expressed by a function of Xj^, X^ and Bq. 
Since the Hamiltonian is expressed by the pair operators X^, X^ and Bq, we can 
rewritten it by the intrinsic pair operators and AN . 

H = {H) + Hrot + Hintr + Hcoupl- (24) 
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Here Hrot is a function of AN 

Hrot = Erot + CiAN + C2AN^ + ■■■, (25a) 
Erot = -C2{AN^) (25b) 

Hintr is composed of only the intrinsic pair operator, 

Hintr = i^inL + ^-f^intr + ^^-^intr ' (^6) 

the leading term of which is given as 

nil = E E-^-a + E Vx{nn')XiXn' + E Mnn') (xt^, + h.c.) , (27) 

<^ nn' nn' 

Vx{nn') = E ■^aP;j5ll^n{aP)lpn'{'yS) 

ajd-fS 

- E(^a + Ep)MaP)^l;r^'{aP), (28a) 

a/3 

Fv(nn') = E ^a/3;75V'n(«/3)V'n'(7<5), (28b) 
a/37<5 

where Ea denotes the quasi-particle energy. H^oupl is a function of AN and intrinsic 
pair operators. It should be noticed that and Xc terms in Eqs.(22) and (23) 
are absorbed in AN term in Eq.(24), because the relation [H,AN] = should be 
hold. Thus we separate the original Hamiltonian into three parts, the rotational 
part, intrinsic part and their coupling part. 

Now we span the intrinsic subspace by a set of intrinsic states 

|0), Xt|0), Xtxt,|0),..., (29) 

which are named generically as 

\k) = Xlxl---Xl\Q). (30) 

Then the number projected intrinsic states are given by 

\k))=PN\k), Pn^ / ^e^"^^^. (31) 

Jo 

The energy matrix and the norm matrix between these states are given by 

Hkk' = {{k\H\k')) 

= (0| • • • Xr^.PNH PnXI^ • • • ^I; |0 ) 

= {0\PN\0)(k\ Hintr \k'), (32) 
Mkk' = fk\k')) 

= {o\Xnr--Xn,PNxl',---K', |0) 

= {0\PN\0){k\k'). (33) 
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Here we use the relation PnAN = AN Pn = and [Pn,XI] 
Therefore the eigenvalue problem reduces to that of Hintr, ic, 



k' 



[PN,Xn] = 0. 

(34) 



where Hkk' = (^| ^mtr Nkk' = {k\k'). Note that (0| P/v |0) is factored out so 
that we need no integration. 

Since the lowest order of Hintr is given by Eq.(27) we can diagonalized it by 
introducing the RPA mode in the intrinsic subspace (which we call IRPA mode) 



5 = E['^A(n)Xt-3;;,(n)X, 



(35) 



(36) 



The wave functions, X\{n), y\{n), are determined so as to satisfy 

[^>a]=^Aa+o(£), 

which gives the IRPA equation for X\{n) and yx{n)\ 

uJxXxin) = EnMn) + ^ Vx{nn')Xx{n') + ^ Vv{nn')yx{n'), (37a) 

n' n' 

-ooxyxin) = EMn) + ^ Vx{nn')yx{n') + ^ Vv{nn')Xx{n'). (37b) 



Thus H^ntr written as 



A 

Eit = -E^.yl{n). 



(38) 
(39) 



It should be noticed that the IRPA vacuum is defined as the vacuum for only the 
intrinsic mode and no difficulty arises in the conventional RPA vacuum which is 
defined as the vacuum for all the eigen mode including the zero-energy mode. 

The higher order terms of Hintr are given by more complicated functions of Xj^, 
Xji, Bq. We must generally solve the eigenvalue problem of Eq.(34) for the higher 
order calculations. 

In order to illustrate the appli- 
cability of our method we calculate 
the ground state energy of the sys- 
tem with the constant pairing in- 
teraction. It is well known that the 
exact solution for the state with se- 
niority is easily obtained numer- 
ically. The numerical results for 
Sn isotopes are shown in Table I. 
The single particle energy and 
the pairing force strength are taken 



Table I. Ground state energies (in MeV) for Sn iso- 
topes are shown. 



A 


BCS 


RPA 


PWO 


PW2 


exact 


104 


-2.313 


-2.742 


-2.742 


-2.744 


-2.745 


108 


-1.844 


-2.526 


-2.526 


-2.534 


-2.534 


112 


0.930 


0.140 


0.140 


0.131 


0.134 


116 


6.030 


5.216 


5.216 


5.193 


5.180 


120 


13.489 


12.761 


12.761 


12.758 


12.761 


124 


22.903 


22.330 


22.330 


22.334 


22.333 


128 


33.837 


33.499 


33.499 


33.501 


33.500 
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from Ref. We performed the calculatfon up to O(e^) and refer it as PW2, while 
the calculation up to 0(e°) is denoted by PWO. We find that the results of PWO 
arc completely agree with those of the conventional RPA. This is not trivial because 
each value of Erot and Eyih is different between RPA and PWO. It is also found that 
PW2 calculation generally improve the RPA results. 

In summary we propose a new method to treat the number conservation. By in- 
troducing the intrinsic pair modes which commute with the particle number operator 
the original Hamiltonian is divided into three parts; the rotational, the intrinsic and 
the coupling part. Thereby we never treat the rotational mode as the zero energy 
eigen mode, and therefore we circumvent the difficulties lying in the conventional 
RPA. The wave function is formally projected to the states with definite particle 
number. However we need no integration over the phase angle as in the Tisual pro- 
jection method because the eigen value problem is reduced to that of the intrinsic 
subspace. The numerical calculation was done to illustrate the applicability of our 
method and the result shows that our method works well. The detail will appear in 
a subsequent paper. 
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